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S£ Abstract 

In a dense cloud of massive fermions interacting by exchange of a light scalar field, the effective 

O 

mass of the fermion can become negligibly small. As the cloud expands, the effective mass and 
j> the total energy density eventually increase with decreasing density. In this regime, the pressure- 

density relation can approximate that required for dark energy. We apply this phenomenon to the 
expansion of the Universe with a very light scalar field and infer relations between the parame- 
ters available and cosmological observations. Majorana neutrinos at a mass that may have been 
recently determined, and fermions such as the Lightest Supersymmetric Particle (LSP) may both 
be consistent with current observations of dark energy. 



tgoldman@lanl.gov 



gjs@phys.unm.edu 



* Electronic address: 
^Electronic address: 
^Electronic address: 
^Electronic address: bhjm@unimelb.edu.au 



alsingpm@hotmail.com 



I. INTRODUCTION 



Several years ago it was suggested that neutrinos might interact weakly among themselves 
through the exchange of a very light scalar particle [TJ [2], with possible consequences for 
the evolution of the Universe and for the propagation of neutrinos from distant events. We 
examined such a system for scalars with astrophysical ranges to explore the possibility of 
neutrino clustering [3J and noted at the time that the neutrino clouds thus formed could 
seed structure formation in the Early Universe. More generally, in such clouds of massive 
fermions interacting by exchange of a light scalar field, the effective mass of the fermion can 
become negligibly small. We found that, as a consequence, when the cloud expands, the 
effective mass and the total energy density must eventually increase with decreasing density. 
We studied this system in 1996 [3] , well before the discovery of Dark Energy, in connection 
with experimental problems encountered in the search for the mass of the (electron) neutrino. 
Those anomalies have since disappeared, but provided us with the technology to describe 
dark energy in a well understood dynamical system. 

In the following, we first review our previous work on the theory of massive fermions 
interacting via exchange of a scalar field. This is carried out with scaled variables so the 
regime of applicability is not constrained. We next recall the relation between Dark Energy 
and equations of state and define the w parameter used therein. After this, we apply our 
model results for w and discuss the numerical, analytical and scaling properties relevant to 
the accuracy of our results. Penultimately, we extract a rough mass value from applying 
those results to describe Dark Energy assuming the currently accepted value for its energy 
density in the epoch corresponding to z — 1. Finally, we present our conclusions and discuss 
some open questions. 

II. SUMMARY OF A THEORY OF MASSIVE FERMIONS 
INTERACTING VIA LIGHT SCALAR FIELD EXCHANGE 

The effective Lagrangian for a Dirac field, ip, interacting with a scalar field, 0, is: 
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which gives as the equations of motion 
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As usual, we set h — c — 1. We have omitted nonlinear scalar self couplings here, even though 
they are required to exist by field theoretic selfconsistency, jl] as they may consistently be 
assumed to be sufficiently weak as to be totally irrelevant. The parameter rrif is the 
renormalized vacuum mass that the fermion would have in isolation, and takes into account 
any contributions from all other interactions, as well as contributions from the vacuum 
expectation value of the new scalar field, <fi. 

We look for solutions of these equations in infinite matter which are static and transla- 
tionally invariant. Eq.([2]) then gives 
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which, when substituted in Eq.([3]) gives an effective mass for the fermion of 
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These equations are operator equations. We next act with each of these equations on a 
state \ defined as a filled Fermi sea, with a number density n per fermion state, and Fermi 
momentum fcp, related as usual by n = /cf 3 /(67t 2 ). The operator ipip acting on this state 
gives 
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where ( is the number of fermion states which contribute — £ = 2 for Majorana fermions 
and ( = 4 for Dirac fermions. Thus the effective mass is determined by an integral equation 

9 2 C ^ 



m 



f 



m 



(0) 

/ 



2ir 2 m 2 Jo 



k z dk 



ni 



f 



m 



f) 



+ k 2 



(7) 



To discuss the solutions of this equation, we reduce it to dimensionless form, dividing by 
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and introducing the parameter 



(OK 2 



.9 2 {rnf) 



2tt 2 



2 ' 



(8) 



and the variables y 
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Then Eq.(7) becomes 
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with ep = yxp + y 2 . This choice of scaled variables gives all energies (and momenta) in 
units of the vacuum fermion mass. For consistency, we define the dimensionless scalar mass 

as a non-linear equation for y as a 
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The total energy of the system is a sum of the energy of the fermions, Ef — e/ (N, 



and the energy in the scalar field, E s 
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(N, where N is the total number of neutrinos 



in each contributing state. These expressions serve to define the per fermion quantities ef 
and e s . Also, E s = £ s V, where £ s = hm 2 (j) 2 is the energy density of the (here uniform) 



scalar field. 

One finds that 



xp Jo 



x dxJx 2 + y 2 



X' 



3_ \ e F xp e F x F y 2 _ j^. 
4 ' 8 8 11 



cf + xf 

. y , 



(ii) 



and 



2 xf 
1 3 



X F 



1 



o \/x 2 + y 2 1 

yf- 



and the total energy density per fermion is just the sum, 

< e >= e/ + e s . 

Notice that for large values of xf, 
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It is also useful to note that, for small xf, 
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For the fermion system to be bound, the minimum of e = Cf + e s as a function of 
density (or :rp) must be less than 1, its value in the zero density limit. Fig.Q shows the 
variation of e and y as a function of xp for several values of Kq. Note that for sufficiently 
large Kq, there is a minimum relative to both the large and small xf regimes, that is, relative 
to regions of both large and small fermion density. 

Thinking of this in terms of expansion of the Universe, early times correspond to the 
large xf region on the right and late times, including presumably the present, are to the 
left. Thus we see that at some intermediate period, the energy density passes through a 
minimum and as the system approaches the present, it passes through a regime in which 
the energy density is increasing as the number density decreases - characteristic of a regime 
of negative pressure. Fig.Q gives an advance peek at the value of the equation of state 
parameter, w, that we derive from the dependence of e vs. xp- We will demonstrate later 
how we do this numerically, but it should be noted that there are still some numerical 
difficulties evidenced by the "hash" in the low xf limit where we know that w — > as the 
cold, now non-interacting fermion "dust" turns effectively into new "dark matter". 



III. EINSTEIN, FRW AND EQUATIONS OF STATE 

In a Friedmann-LeMaitre-Robertson- Walker Universe, Einstein's equations produce the 
second order time derivatiive equation of motion that relates the expansion (size) scale 
parameter, a, to Newton's constant, G, the matter density, p, pressure, P, and a cosmological 
constant, A: 

a 4:7rG „ s A , . 

- = ~j- (p + 3P)+3 (16) 

Therefore, it is necessary to know the relevant equation of state (EoS) before the time 
development of the scale factor can be determined. For dust, which has no pressure, P = 0, 
while for a relativistic gas, P = |p. Note that acceleration of the expansion parameter 
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FIG. 1: Total energy density per fermion e and effective mass y vs. xf for several values of Kq. 

occurs, even in the absence of a cosmological constant, when P < — |p. Finally, for a 
spatially and temporally homogenous scalar field, P < —p, which is more than enough to 
produce acceleration of the expansion. More generally, we can parametrize equations of 
state in this regime by a constant, w, as 

P = wp (17) 

There are a great number of models for this Dark Energy phenomenon. They go by 
names such as "quintessence" for — 1 < w < 0, or "phantom energy" for w < —1. Among 
others, this issue has been addressed by Fardon, Nelson and Weiner [6], Peccei [7], Barshay 
and Kreyerhoff [8], Baushev [9] and Mukhopadhyay, Ray and Choudhoury [10J. 



A. Our EoS 



For the system under consideration, the total (matter plus field) energy is given by the 
product of the total energy density per fermion < e > and the total number of fermions, 
which in turn is determined by the number density n times the volume, V: 



U = P V = mf ) <e> nV, 



where we recall that 
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FIG. 2: y and w vs. log(xp) for Kq = 10 
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where U is the internal energy given above. Since nV is a constant, 
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From this and Eq.(17), we can identify 
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This is our central model result. 



IV. GENERAL CHARACTER OF MODEL RESULTS 

In Fig.Q, we shows the value of w as computed numerically from Eq.(22) for 8 values 
of K on a log scale for xp < 0.4 and in FIg.Q for < 0.1 on a linear scale. Note that 
it approaches close to —1 as the density decreases (as the Universe expands and the scale 
factor a increases from right to left) and then departs sharply towards zero, as also indicated 
earlier in Fig. ([2]). 

At large xf, it is clear that w approaches +1/3 as it should for a relativistic gas of 
fermions. It is perhaps less clear, due to numerical fluctuations, that the value goes to zero 
at zero density. This can be checked by considering the small xp expansions of the energy 



densities shown at Eq.(15). 



Finally, we performed a number of numerical scaling checks to examine whether w can 
fall below —1. The are shown in Figs.( 5p|7 ). 



We are continuing our efforts to ensure numerical stability primarily by seeking analytic 
formulae and approximations, especially to avoid taking derivatives numerically. We will 
report on these improvements elsewhere, but suffice it to say that we have found no con- 
tradiction to the results obtained here by purely numerical means, and confirmed that the 
difficulties encountered at very small xp are indeed due to numerical noise. 
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FIG. 3: w vs. log(x^) for 8 values of Kq. 
V. NUMERICAL VALUES 

All of the above is carried out with scaled variables. However, was shown in Ref. [3], there 
are only a few, weak constraints on the actual parameter values. Very large values of Kq are 
possible even for very small values of g 2 if the range of the scalar is very large, corresponding 
to very small values of m s . Even if long-ranged, such weak interactions between fermions, 
especially neutrinos or those outside of the Standard Model altogether (such as the LSP) 
are exceptionally difficult to constrain by any laboratory experiments. 

The energy density of Dark Energ y is quoted [5J as (3.20 ± 0.4) x 10~ 47 Gev 4 . In more 
manageable units, this is given as (2.3meV) 4: . What does this imply for the allowed value 
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FIG. 4: w vs. xf for 8 values of Kq. 

of mo? If we set this energy density equal to that of this system at w = — 1, 

(m ( f 0) ) 4 < e > xi 7 , 



then solving for mj ' gives 



If we further suppose that this occurs at cosmo logical z ~ 1, then the range of the scalar 
field must be comparable to the size of the Universe at that time. That is, 

m s ~7x 10 9 lightyears ~3x KT 30 meV. (25) 

For the relatively modest value of K ~ 8 x 10 6 , this implies that mf^ ~ 300meV and 
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FIG. 5: Deviation of minimum value of w from —1 vs. Kq. 

g 2 /(47r) ~ 6 x 1CT 58 . These values emphasize the virtual impossibility of constraining this 
physics by means of laboratory experiments. 

We note with interest that following the curve for y from z — 1 to "now", i.e., z = 0, 
tells us that the effective mass of this fermion would now be measured to be approximately 
(7/8)mj°' ) , a value tantalizingly close to the Majorana neutrino mass that Prof. Klapdor has 
reported from his experiments [TT] . 

Other solutions are possible, and Fig. ^ shows how our results scale very accurately (for 
sufficiently large Kq) with the 4th root of K , both numerically and under one of our analytic 
approximations to the region where w is a minimum. In particular, if one has the LSP at 
~ 1 TeV in mind, Kq becomes very large, ~ 10 57 and g 2 /(47r) also increases, ~ 10~ 32 , but 




FIG. 6: Minimum value of w vs. xf at which minimum occurs, for several values of Kq. 
these values are not ruled out by anything known. 

VI. CONCLUSIONS AND QUESTIONS 

We have displayed an explicit and calculable dynamical mechanism that describes Dark 
Energy and connects it to what in the current epoch becomes a kind of dark matter. Al- 
though Dark Matter is usually thought of as existing at early epochs in the life of the 
Universe, the system described here turns (hot) relativistic fermions, that interact weakly 
with a very light scalar field, into Dark Energy which lasts for a limited time during the 
expansion of the Universe which then morphs into new cold dark matter components. Thus, 




FIG. 7: Minimum value of w vs. logarithm of xp at which minimum occurs, for several values of 
K . 

neutrinos can contribute to both. Nor need there be only one time scale or one species for 
which this applies. (See, e.g., Ref.[12].) If there are many sterile fermions with sufficiently 
long decay lifetimes, the acceleration/deceleration history of the Universe could be much 
more complicated that presently envisioned. 

We may also ask generally why w > — 1, but it is fairly clear in this model: A scalar field 
strength uniform in space and time produces [13] exactly w — — 1, but here the source for the 
scalar field is the density of massive fermions. Relativistically, their strength for producing 
scalar field is severely reduced at high momentum (in the rest frame of the Universe) but 
as they slow, the scalar field strength grows nonlinearly until, due to the expansion of the 
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FIG. 8: Scaling variation of the vacuum fermion mass, , vs. Kq. 

Universe, the fermions separate so much (greater than the Yukawa range for exchange of the 
scalar field) that they cannot act collectively and the scalar field strength declines again. 

Finally, we note that our model has definitive if difficult tests, as it predicts specific 
variations of w, slowly approaching —1 from above as z decreases through 1, and rising 
rapidly towards zero as z approaches the present. We hope this encourages observationalists 
in their efforts to discern variation of w with z. 
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